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1> Abstract 
43 ■ 

We extend the MRS (A) set of parton distributions, which provides an up- 
■ to-date description of the structure of the proton probed by hard interactions, 

into the region of Q 2 below 5 GeV 2 . After physically-motivated modifications, 
we obtain a smooth description of the data on the structure function F2(x, Q ) 
all the way from Q 2 = 10" 1 to 10 3 GeV 2 . 



1 Introduction 



Perturbative QCD is remarkably successful in describing a broad sweep of hard- 
scattering processes involving the nucleon. This implies that it is possible to extract 
from the data a consistent set of parton distributions in the proton, fi(x, Q 2 ), which 
evolve in Q 2 according to the standard GLAP equations in the region where pertur- 
bative QCD is appropriate. 

Experimentally there has been steady progress, first in extending the kinematic 
range and in improving the precision of the data (in particular of the deep-inelastic 
structure function measurements) and, secondly, in extending the range of "hard" 
scattering data (e.g. the recent measurements of the W ± rapidity asymmetry |J 
and of the asymmetry in Drell-Yan production in pp and pn collisions [0]). As a 
result, modern sets of parton distributions || |], |5| are strongly constrained. Here we 
concentrate on the most recent set of partons, MRS (A) ||, which were obtained in a 
global analysis of the full range of deep-inelastic and related data. These distributions 
provide a detailed description of the proton structure in the region where leading- 
twist perturbative QCD is valid, which may be taken as Q 2 > Ql ~ 4 GeV 2 . In 
the case of deep-inelastic scattering, we know that higher-twist contributions become 
important even above this value of Q 2 , especially at large x, and so an additional cut 
is always made in W 2 , typically W 2 > W$ ~ 10 GeV 2 . 

There are, however, good reasons for requiring a quantitative description of the 
nucleon structure at lower values of Q 2 . In structure function measurements, for 
example, the implementation of radiative corrections requires an approximate de- 
scription of F 2 at low Q 2 . Indeed, there now exist several low-Q 2 parametric forms 
for F 2 (x, Q 2 ) |j, 0, H, |9], [nj which accurately describe F 2 in a restricted region of 
Q 2 and x. An excellent summary of the situation has recently been presented by 
Badelek and Kwiecinski JTI]. There is also a demand for individual parton distri- 
butions fi(x,Q 2 ) valid at low Q 2 as well as for Q 2 > Qq. For example, many of 
the measurements at HERA involve the use of Monte Carlo event generators which 
describe parton showering down to really low Q 2 values. 

The aim of this paper is to provide a set of parton distributions that is consistent 
with the data taken at low Q 2 (down to Q 2 ~ 0.1 GeV 2 ) and, as Q 2 rises, smoothly 
approaches the set MRS (A) which is consistent with the high-Q 2 experimental data. 
Note that there is already one set of parton distributions, that of Gluck, Reya and 
Vogt (GRV) fli"2| , which could, in principle, be used down to Q 2 = 0.3 GeV 2 . Their 
procedure is to start with a set of 'valence'-like quark and gluon distributions at fi 2 = 
0.3 GeV 2 and to evolve these up in Q 2 using the standard leading-twist (= 2) next-to- 
leading-order evolution equations. Above Q 2 = 0(2 GeV 2 ) these provide a reasonable 
description of F 2 , consistent with the belief that, by then, the twist-two component 
is expected to dominate. At low Q 2 , however, the higher-twist contributions are 
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important and the GRV partons alone do not give a good description of F 2 there, as 
was never their intention. Moreover, there are also concerns that the perturbative 
GLAP evolution of valence-like GRV partons at small x and Q 2 is unreliable |T3] . 



In contrast, what we provide here are 'effective' low-Q 2 parton distributions which 
are modified versions of the leading-twist distributions, i.e. they satisfy leading-twist 
Q 2 evolution only for Q 2 > 5 GeV 2 . In other words, the unknown higher-twist and 
non-perturbative components are the origin of the empirical Q 2 modifications we im- 
pose on the distributions at low Q 2 , such that when these partons are inserted into 
the formal expressions for F 2 , we obtain detailed agreement with the data. Particu- 
larly useful in this context are the new low-Q 2 structure function measurements from 



the E665 collaboration 14 



Note that the gluon and sea-quark distributions of the MRS (A) set of partons 
have the singular small-x behaviour 



xg, xq sea ~ x ■ (1) 

as x — > 0, which results from including the HERA measurements of F 2 in the global 
analysis. For such a singular gluon distribution, GLAP evolution is expected [13|, [15 
to faithfully mimic the small-x description of F 2 by the BFKL equation [TH] . 



2 Parametrization of partons at low Q 2 

Our starting point is the MRS (A) set of parton distributions, which are obtained 
from a global data analysis of deep-inelastic and related data with Q 2 > 5 GeV 2 (and 
W 2 > 10 GeV 2 ). The partons give an excellent description of these data. However 
our goal here is to obtain partons which describe data to much lower values of Q 2 
(and W 2 ). We therefore begin by using the next-to-leading-order GLAP equations 
to evolve down in Q 2 from the MRS (A) starting distributions at Ql = 4 GeV 2 . At 
some stage we anticipate that such a leading-twist GLAP-based extrapolation will 
fail. Note that as we go below the charm threshold, the number of active flavours in 
the evolution drops from 4 to 3, and from our previous analysis || this threshold is 
taken at Q 2 = m 2 c = 2.7 GeV 2 . 

First we must allow for the effects of the target mass, which were not included 
in the MRS (A) analysis. Such effects do not so much concern low-x data (where 
the photon-proton scattering energy is large) but rather apply to data at large x 
where they are significant even for Q 2 > 10 GeV 2 . Data in this kinematic region 
were excluded from the MRS (A) analysis by the W 2 cut and so did not distort the 
fit. To include the target mass corrections we replace x by the target mass variable 
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2x .., 2 4M¥ 



£ |T7|, where 

with r> = l + ^f- (2) 

where M is the mass of the target. From Fig. 1 we see that this straightforward 
correction gives a dramatic improvement in the description of the large x SLAC 
measurements of F 2 ||18|| . Here the difference between the modified and unmodified 
data is almost entirely due to the x — > £ substitution. 

Our primary interest is the behaviour at low Q 2 . There exist fixed-target mea- 
surements of F 2 (x,Q 2 ) at low Q 2 and low x. Indeed, the E665 collaboration have 
recently presented preliminary data which access even lower Q 2 and x values than 
hitherto. We expect that F 2 (x, Q 2 ) reconstructed from MRS(A) partons evolved to 
low Q 2 will, at some stage in the backwards evolution, begin to overshoot the data, 
since the fundamental requirement that 

F 2 - 0(Q 2 ) (3) 

as Q 2 — > is not embodied in the perturbative QCD evolution. The most natural 
way to rectify this deficiency is to introduce a form factor 

QZ _|_ m Z 

into the parton distributions, where the x dependence of m 2 is to be determined from 
the low-Q 2 data. In this way we achieve a phenomenological description of higher- 
twist contributions which may, for example, arise from parton shadowing [20]. If 



the latter is the dominant higher- twist effect, we would expect m 2 to increase with 
decreasing x, as can be seen from the following simple argument. Naively, we would 
expect shadowing to become important when the total gg 'interaction area' becomes 
a significant fraction of the transverse area (vri? 2 ) in which the gluons are confined 
within the proton. That is when 



ng(J " ~ 0(a s ) , (5) 



7TR 2 

where n g ~ xg is the gluon density per unit rapidity. Noting that a gg ~ a 2 /Q 2 , this 
implies that the 'critical line' in the (Q 2 , 1/x) plot for the onset of shadowing has the 
form 

Q 2 ~ m 2 „ x - x (6) 

where, for simplicity, we have assumed xg,xq ~ x~ x and ignored the running of a s . 

In practice it is clearly desirable to distort the MRS (A) partons as little as possible 
for Q 2 > Qq = 4 GeV 2 , and so we choose 

m 2 = m 2 (x) exp(-Q 2 /Q 2 ) . (7) 
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Thus there is negligible modification to the partons (i.e. p ~ 1) in the kinematic 
region (Q 2 > Ql) of the data that originally determined the MRS (A) distributions. 
The exponential factor in (0) does not significantly distort the determination of the 
x dependence of m 2 , since m 2 is mainly determined by F 2 data with Q 2 < 1 GeV 2 . 

In summary, we modify the evolved MRS (A) parton distributions so that they 
take the form 

f t (x,Q 2 )=p(x,Q 2 )fr SA (^Q 2 ), (8) 

and determine the parameter wig in p by fitting to the F 2 (x, Q 2 ) data at low Q 2 for 
each value of x. It is interesting to note that the form factor modification and the 
x —>■ £ target mass correction essentially decouple as they are relevant in distinct 
kinematic regimes, namely in low-Q 2 , low- a; and in large- a; , moderate-Q 2 regions 
respectively. 

In principle it could be argued that a different form factor should be used for 
different parton distributions. However the low-Q 2 , low-x data are dominantly de- 
scribed by the sea-quark distributions. For instance, we find no significant change in 
the results if we use an unmodified distribution for the gluon, that is p = 1. (Com- 
puting the value of R = o L jo T using the QCD expression for F L evaluated with our 
modified distributions does not remove the discrepancy between the SLAC data 



and the QCD prediction. The values of R obtained by SLAC still tend to lie above 
the computed values.) 

Finally we note a small technical point. When the MRS (A) parton distributions 
are evolved down from Ql = 4 GeV 2 , the gluon distribution xg(x,Q 2 ) begins to 
become negative over a small range of moderate x, just below Q 2 = 0.625 GeV 2 . We 
therefore freeze the unmodified partons, f™ RSA (x,Q 2 ), for Q 2 < 0.625 GeV 2 at their 
values at Q 2 = 0.625 GeV 2 , i.e. 

Q 2 - Q 2 + (Ql - Q 2 )0{Q 2 c - Q 2 ) (9) 

in the argument of the ff AnSA , with Q 2 = 0.625 GeV 2 . The effect on the behaviour 
of the u sea-quark distribution is shown in Fig. 2 for two relevant values of x. Since 
the unmodified distributions are decreasing only slowly with Q 2 there is little effect 
on the description of F 2 . 



3 The description of data at low Q 2 

The behaviour of the parameter m 2 (x) as a function of x is determined by fitting 
to the F 2 data of E665 and SLAC [0, and also the 90 GeV data of NMC [[ZTJ. 
The NMC data for Q 2 > 5 GeV 2 proved to be crucial in constraining the MRS (A) 
partons. Although in principle the overall normalization of the NMC data could be 
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varied with respect to the other data sets, no correction was in fact needed ||. When 
we include the E665 data in our analysis we again allow for an arbitrary relative 
normalization adjustment. There is in fact a significant overlap in the x, Q 2 ranges 
of the NMC and E665 data sets and we find that consistency is achieved if the latter 
are renormalized up by 20%. Although this may appear to be a large factor, the 
E665 data are still preliminary, with a figure of 10-20% being quoted as the typical 
systematic uncertainty [14|. Applying this correction factor 1.2 to the E665 data leads 



to a systematic x dependence in the values of m^x) extracted from F 2 measurements 
of the three experiments, SLAC, NMC and E665. 

The results are shown on a logm 2 , — log a: plot in Fig. 3(a), and indicate that a 
good fit to the x dependence of ml can be obtained by using the form 



m {x) 



Ax~ 



(10) 



We find A = 0.07 GeV 2 and n = 0.37, corresponding to the straight line in Fig. 3(a). 
It is remarkable that the x~ - 37 behaviour obtained from low-Q 2 data is compatible 
with the observed singular x -0 3 gluon and sea-quark small- a: behaviour, flip, obtained 
from large Q 2 data. This connection between independent results obtained from 
different kinematic regions is suggested by the simple arguments which lead to (M). 



It is interesting to note that 
alternative form 



(x) could be equally well by represented by the 



rrir 



(x) = Aexp 




(11) 



as shown in Fig. 3(b). The straight line corresponds to A = 0.015 GeV 2 and B = 1.54 
in (|TT|). Note that the quoted values for A and n in ( [TD[ ) are for fits are performed in 
the MS scheme. A similar fit in the DIS factorization scheme gives A = 0.055 GeV 2 
and n = 0.39. 

In Figs. 4 and 5 we show the description of the low-Q 2 measurements of F 2 |b|, [18 



21 1 using the MRS (A) partons modified as in (||) with m\ given by (]T0f) . In Fig. 6 we 
show, at sample values of x, the continuation of the low-Q 2 description to the recent 
HERA measurements of F 2 [22, [23|. The continuous curves show that the partons 
give a good description of F 2 throughout the range 10 _1 < Q 2 < 10 3 GeV 2 . The 



predictions of the 'dynamical' GRV partons [|12[] are also shown (by dotted curves) in 
Fig. 6. Within the scope of their model, these partons provide an excellent description 
of F 2 for Q 2 > 1 GeV 2 but below this value increasingly undershoot the data. It can 
also be glimpsed from Fig. 6 that the recent HERA data indicate that the GRV 
predictions increase a little too steeply with decreasing x. 
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4 F 2 n at low Q 2 



There is also evidence for higher-twist corrections from SLAC and NMC low-Q 2 data 
24 , |25[| for the ratio of the neutron and proton structure functions. The deviations 



from the MRS(A) perturbative QCD predictions for F-f/Ff occur in the interval 
x ~ 0.1 to 0.3 and so we introduce a further modification to the valence quarks solely 
to describe these deviations.^ We use the following simple one-parameter form which 
allows F£ to be varied at low Q 2 while leaving Ff unaffected: 

u v (x, Q 2 ) -> u' v (x, Q 2 ) = (1 - ^)u v (x, Q 2 ) - ^d v (x, Q 2 ) 

d v {x,Q 2 ) -> d' v (x, Q 2 ) = (1 + r)d v (x, Q 2 ) + ru v (x, Q 2 ) (12) 

where 

r = Q \ . (13) 

The F^/Ff data obtained by NMC are shown in Fig. 7. They have been corrected 
for deuteron shadowing effects using the results of the analysis of |f26[. We see that the 
predictions obtained from the MRS (A) partons (continuous curves) tend to under- 
shoot the data. This discrepancy, however, is much less than for Ff itself at low Q 2 , as 
would be expected for the ratio F^ / Fj which depends primarily on the valence distri- 
butions. The dashed curves, obtained using Eqs. ( |T2"D and flTB| ) with Q 2 = 0.12 GeV 2 , 



give a good description of the data. Whether this simple parametrization is adequate 
for smaller values of Q 2 is not certain. In any case, nuclear shadowing at such low 
Q 2 is expected to overwhelm such a simple partonic description. 

In Fig. 8 we show the prediction of the modified MRS (A) partons (continuous 
curves) for the EMC(NA28) [f^J] measurement of F 2 on a deuterium target. In this 
case the curves (and not the data) have been corrected for deuterium screening effects. 
The description is very satisfactory and mainly checks the form factor modification 
of (|]). These data are not sufficiently precise to further constrain the parameter r in 
(|T^), since the valence-quark contribution is small for those EMC data which lie at low 
Q 2 . For comparison, we also show the GRV predictions for F 2 D . At Q 2 = 0.35 GeV 2 
the GRV values are about a factor 2 below the data, but very rapidly evolve upwards 
with increasing Q 2 to be in reasonable agreement with the data. 



1 Note that the flavour-independent form factor modification of Eq. (^) would leave the F-fV-Ff 
ratio essentially unchanged. 
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5 Conclusions 



A low-Q 2 parton model is really a contradiction in terms. At some stage as Q 2 de- 
creases, the description of the proton's structure cannot be expressed in terms of 
single parton densities with simple logarithmic behaviour in Q 2 . The contributions 
from parton-parton correlation densities with power behaviour must enter and even- 
tually the non-perturbative behaviour dominates. Moreover, all of these contributions 
are expected to be process-dependent. The set of distributions we have derived in 
this study correspond to what we may call 'effective' partons, insofar as power-law 
corrections like those introduced in Eq. @ would presumably arise from including 
gluon-quark or quark- ant iquark correlation distribution functions. In the framework 
of the operator product expansion, we assume that these higher-twist terms can be 
added to the leading-twist piece, each with their separate Q 2 dependence. These 
effective partons do not of course satisfy the usual conservation laws - number of 
valence quarks, total fractional momentum, etc. - but the violation is below 20% 
even for Q 2 ~ 0.5 GeV 2 . Their detailed structure at low Q 2 is also very dependent 
on the precise form of the starting distributions at Qq = 4 GeV 2 . Note that, in 
contrast with the studies in Refs. |], [7|, |S], [|, [TIJ , we are not attempting here to make 
a model which has a smooth transition to the Q 2 — > photoproduction limit, where 
a parton-based approach is clearly invalid. 

From our study we may estimate, at least in the case of 'singular' MRS-type 
parametrizations, where the leading-twist GLAP evolution begins to become unreli- 
able as Q 2 decreases. Interestingly, this depends on the value of x in just such a way 
as if the higher-twist contribution arose from parton shadowing. Some insight may be 
obtained by presenting the results as a function of x at fixed Q 2 . In Fig. 9 we display 
this dependence of the F% and _F 2 D data for Q 2 ~ 0.3 GeV 2 . In this low-Q 2 regime 
the data are reasonably flat, whereas we see that the unmodified MRS (A) partons 
(dashed curves) partly retain the x -0 ' 3 behaviour of the distributions at the start 
of the backward evolution. The higher- twist, form-factor modification (continuous 
curves) restores the agreement with the data, with the modification decreasing as x 
increases. These data indicate large higher-twist contributions, but as Q 2 increases 
towards 5 GeV 2 these effects rapidly decrease. 

Backward evolution in Q 2 is much more sensitive to the starting distributions 
than is forward evolution. It is therefore possible to conceive 'non-singular' starting 
parametrizations (with a lower starting scale Qo) in which the higher-twist effects 
estimated in this way are much smaller. In this case the rise in F 2 with decreasing x 
seen at HERA is generated by the long evolution length. The GRV partons |12| are 
a parametrization of this type. We show in Fig. 9 the GRV predictions, which at this 
value of Q 2 undershoot the data with a shape which reflects their valence-like input 
at Q 2 = 0.3 GeV 2 . With increasing Q 2 , the agreement between the GRV description 
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and the data rapidly improves (see Fig. 6). 

Our justification for discussing individual parton distributions at very low Q 2 
values is largely practical. As mentioned earlier, parton showering Monte-Carlo pro- 
grams at HERA incorporate parton densities at some low cut-off in Q 2 . We are 
providing here effective parton densities which, when substituted into the next-to- 
leading-order expressions for F 2 (x,Q 2 ), give a reliable description of the data down 
to Q 2 = KT 1 GeV 2 and, on the other hand, up to Q 2 = 10 3 GeV 2 .[] Fig. 6 shows 
how F 2 (x, Q 2 ) for small values of x is well described throughout all of this region. 
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Figure Captions 



[1] Measurements of F%(x, Q 2 ) at large x obtained by SLAC @ and BCDMS [[19 
The latter data are normalized by a factor 0.98 as required by the global parton 
analysis. The broken curves are obtained using MRS (A) partons H, and the 
continuous curves using MRS (A) partons modified as in (f|). For these data the 
target mass correction (the substitution x — > £) is the dominant modification. 

[2] The sea-quark distribution xu(x, Q 2 ) at two values of x. The upper curves cor- 
respond to MRS(A) partons, frozen at Q 2 = 0.625 GeV 2 . The lower continuous 
curves correspond to MRS (A) partons with the form factor modification as in 
(f|) and (H). Also shown for comparison (dashed curves) are the predictions of 
the GRV set of partons [|12j. 



[3] (a) The values of m 2 obtained from fitting to E665 @, SLAC pf and NMC 
[ PHI measurements of F$(x, Q 2 ) at different values of x using MRS (A) partons 
modified as in (|J). The straight line, m\ = 0.07s -0 37 GeV 2 , is the least-squares 
fit to the values of m^. 



(b) As for (a), but with logm 2 , plotted as a function of J\og{l/x). The straight 



line is the least squares fit, m 2 , = 0.015 exp[1.54'Jlog(l/a:)] GeV 2 



[4] The description of the E665 measurements \TA\ of F^(x, Q ) by the MRS(A) 
partons modified as in (|8]). The curves have been renormalized downwards by 
a factor 1.2 for the reasons described in the text. 



[5] The description of the SLAC Jl8[ and NMC [21] measurements of F%(x, Q 2 ) by 



the MRS (A) partons modified as in 



[6] Fixed-target Q, gTJ and HERA |22|, |23[ measurements of F%(x, Q 2 ) at selected 
values of x, as described by MRS (A) partons modified as in (§). The predictions 
of the GRV parton distributions [12[ are also shown (by the dotted curves) for 
comparison. 

[7] The description of the NMC measurements of F2/F2 by MRS (A) partons 
(continuous curves) and by MRS (A) partons modified as in (d) - (0) (broken 
curves) . 



[8] The description of the EMC(NA28) measurements |27| of F 2 D by MRS (A) par- 



tons modified as in (§) - ( |13|) (continuous curves). The predictions of the GRV 
parton distributions [jl2[ at very low x are also shown (by the dashed curves) for 
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comparison. The statistical and systematic errors on the data have been com- 
bined in quadrature. There is an additional overall normalization uncertainty 
of ±7%. 

[9] The description of the E665 F% (renormalized upwards by a factor 1.2) [jl4| and 
EMC(NA28) F 2 D |7| @ measurements at fixed Q 2 ~ 0.3 GeV 2 by MRS (A) 



partons modified as in (g) - ( |T3"D (continuous curves). The unmodified MRS (A) 
predictions (dashed curves) and the predictions of the GRV parton distributions 
(dotted curve) are also shown. 
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